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The hydrogen adsorption sites in MOF5 were determined using neutron powder diffraction along 
with first-principles calculations. The metal-oxide cluster is primarily responsible for the adsorption 
while the organic linker plays only a secondary role. Equally important, at low temperatures and 
high-concentration, H2 molecules form unique interlinked high-symmetry nano-clusters with inter- 
molecular distances as small as 3.0 A and H2-uptake as high as 10-wt%. These results hold the key 
to optimizing MOF materials for hydrogen storage applications and also suggest that MOFs can be 
used as templates to create artificial interlinked hydrogen nano-cages with novel properties. 

PACS numbers: 61.66.-f,61.12.-q,68.43.Fg,68.43.Bc 
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The success of future hydrogen and fuel-cell technolo- 
gies is critically dependent upon the discovery of new ma- 
terials that can store large amounts of hydrogen at am- 
bient conditions P, Metal-organic framework (MOF) 
compounds, which consist of metal-oxide clusters con- 
nected by organic linkers, are a relatively new class of 
nano-porous material that show promise for hydrogen 
storage applications because of their tunable pore size 
and functionality H 0, H US H Yet despite nu- 
merous experimental studies of hydrogen adsorption in 
MOF materials, the nature of the MOF-hydrogen inter- 
action and the manner in which hydrogen molecules are 
adsorbed onto the structure are still unknown. Answers 
to these questions hold the key to optimizing these ma- 
terials for practical hydrogen storage applications. 

Here using the difference-Fourier analysis of neutron 
powder diffraction data along with first-principles total- 
energy calculations, we directly determined the H2 ad- 
sorption sites in M0F5 (the most widely studied MOF 
material, which consists of Zn04 clusters linked by 1,4- 
benzenedicarboxylate (BDC)). Surprisingly, the M0F5 
host lattice has the enough space available to hold many 
hydrogen molecules, up to 10 wt-% at low temperatures. 
The Zn04-cluster is responsible for most of the adsorp- 
tion while the organic linker plays only a secondary role. 

Equally important, we find that at high-concentration 
loadings hydrogen molecules form unique 3-dimensional 
(3D) networks of H2 nano-clusters with intermolecu- 
lar distances of 3.0 A, which is significantly shorter 
than the intermolecular distances of 3.6 A in pure solid 
hydrogen^3i Q| • These findings suggest that MOF ma- 
terials can also be used as templates to create artificial, 
interlinked hydrogen nano-cages. Such materials could 
exhibit very unexpected properties due to confinement ef- 
fects and small intermolecular separation, such as metal- 
lic behavior [121. 

Due to the large incoherent cross section of H2, the 
neutron diffraction data were collected on a deuterated- 
M0F5 sainple, which was synthesized as described in de- 
tail in Refill- Because of the large cubic crystallite sizes 
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FIG. 1: The neutron powder diffraction pattern (A=2.08 
A) of the deuterated-MOF5 host lattice at 3.5 K (dots) plus 
the Rietveld refinement (solid), using space group Fm3m and 
a=25.909 A, and difference plot (bottom). The inset shows 
the high-angle portion of the data along with the fractional 
atomic positions and thermal factors. The refinement was 
characterized by x^=l-195. 



obtained from the synthesis, we ground the sample in a 
helium- filled glovebox prior to neutron powder diffraction 
to eliminate the effects of preferred orientation. Figure 1 
shows the diffraction data from the deuterated-M0F5 
host lattice which was obtained on BT-1 at NIST. The 
pattern was taken over a period of 8 h with the 2.5 gram 
sample loaded into a vanadium cell at 3.5 K. The agree- 
ment between the data and refinement is excellent. We 
note that the neutron diffraction data allowed the de- 
termination of the hydrogen atom positions of the BDC 
linker, which were not previously observed. 

Having characterized the MOFS host lattice, we next 
studied the adsorption of hydrogen in M0F5 as a func- 
tion of D2 concentration per formula unit (i.e. 4Zn) d.^. 
The hydrogen loading was achieved by first filling a well- 
known dosing-volume to a target pressure and then ex- 
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FIG. 2: (a)The target versus refined values for the hydrogen 
loading of M0F5 host lattice, (b) The lattice parameter as a 
function of hydrogen loading. 



posing it to the M0F5 sample at 70K. The sample was 
then cooled down to 30K at which point the pressure 
decreased to a negligible value as the D2 was adsorbed. 
Once the system was equilibrated at 30 K, the sample 
was further cooled down to 3.5 K before the measure- 
ments. The sample was loaded with the following con- 
centrations; nD2=4,8,14, 26, 34, and 46 D2 per molecular 
formula (i.e. per 4Zn)[lJ|. We note that one D2/4Zn cor- 
responds to about 0.255-wt% hydrogen uptake. For the 
nD2=46, the final pressure at 30 K was non-zero, and 
we briefly pumped the system to remove free D2. None 
of the structural refinements of the deuterium loaded- 
samples showed any evidence for solid D2, indicating that 
the deuterium was adsorbed onto the M0F5. This was 
further supported by the total amount of hydrogen ob- 
tained from the refinements as shown in Fig. 2(a). Apart 
from the last point where we had to remove unadsorbed 
deuterium gas, the target and refined values for the to- 
tal amount of deuterium molecules adsorbed in M0F5 
lattice are in very good agreement. It is quite interest- 
ing to note that at cryogenic temperatures, the M0F5 
host lattice actually has enough space to hold up to 10- 
wt% hydrogen. If we can find a way to engineer M0F5- 
H2 interactions in these materials to hold the hydrogen 
molecules in the structure at ambient temperatures, one 
could easily store enough hydrogen for practical applica- 
tions. Figure 2(b) shows that there is a small contrac- 
tion of the lattice upon hydrogen loading. We attribute 
this to a small charge transfer (« O.le) to the hydro- 
gen molecules which results in an attractive Madelung 
Coulomb energy. In the rest of this letter, we discuss 
where and how the hydrogen molecules are packed into 
the M0F5 structure as a function of D2-loading. 

Figure 3 shows the diffraction pattern from M0F5 
which was loaded with 4D2/4Zn. In order to locate the 
hydrogen adsorption sites, we first performed a Rietveld 
structural refinement using the model for the the M0F5 
host structure, ignoring the adsorbed D2 molecules. The 
difference-Fourier scattering-length density based upon 
this model, shown in the inset to Fig. 3, was then used 
to locate the adsorbed D2 molecules. The Fourier plot 
clearly shows where the hydrogen molecules are adsorbed 
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FIG. 3: The neutron powder diffraction pattern (A = 2.08A) 
of the MOF5-4D2 at 3.5 K (dots) plus the Rietveld refinement 
(solid) with space group Fm3m and a=25.9097 A, and differ- 
ence plot (bottom). The inset shows the real-space Fourier- 
difference scattering-length density superimposed with the 
MOF5 structure, indicating the location of cup-sites for the 
first hydrogen adsorption (red-yellow-green region). The re- 
finement was characterized by x^=2.736 and R,i,p=5.54%. 



(red- yellow-green region). Isosurfaces of the three- 
dimensional difference-Fourier scattering-length density 
for a loading of 8D2/4Zn are shown in Fig. 4(a). The 
first adsorption sites (blue) are the positions at the cen- 
ter of the three ZnOa triangular faces, which resemble a 
cup shape and were termed the "Cup-site". There are 
four such sites, forming a tetrahedral cluster (blue region 
in Fig. 4a). Having determined the location of the first- 
adsorption sites, we then further refined the structural 
model, explicitly including the D2 molecules at the first 
adsorption site. The positions, isotropic-thermal factors, 
and fractional occupancies of the adsorbed D2 molecules 
were refined. The deuterium molecules were treated as 
point scatterers with double occupancy. The final re- 
finement for a deuterium loading of 4D2/4Zn is shown 
in Fig. 3. The agreement between data and the fit is 
very good. For the 4D2/4Zn loading, we also observed 
a small amount of D2 (i.e. 10% ) adsorbed at a sec- 
ondary adsorption sites (green isosurface in Fig. 4(a)). 
For 8D2/4Zn loading, these two adsorption sites are al- 
most fully occupied J^]. Unlike the cup-sites, the sec- 
ond adsorption sites are on top of a single ZnOa trian- 
gles and were hence denoted as the "ZnOs-site". These 
sites also form a tetrahedron about the metal-oxide clus- 
ter. We observed that with further hydrogen loading (i.e. 
14 and 26 D2/4Zn), there are two additional adsorption 
sites which start to populate in almost equal proportion. 
These sites are shown in Fig. 4(b) as light-blue and brown 
spheres. The adsorption sites just above the two oxy- 
gen ions are called the "Zn02-site" (see Fig. 4(b)). The 
fourth adsorption site is basically the top of the hexago- 
nal linkers, which we termed the "Hex-site" . The refined 
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FIG. 4: The hydrogen absorption sites obtained from differ- 
ence Fourier analysis. Top: The first (blue) and second (green) 
absorption sites, respectively. Bottom: A view along the 3- 
fold axis, showing the four absorption sites together. In addi- 
tion to the first two absorption sites (blue and green), Zn02 
(brown) and Hex-sites (light-blue) are shown. 



fractional positions of these four sites are summarized in 
Fig. 4. At 26 D2/4Zn loading, structural refinement in- 
dicates that these four adsorption sites are almost totally 
occupied Ts'l, yielding 6.63-wt% deuterium uptake. 

It is important to know if the adsorption sites reported 
above make sense in terms of hydrogen host-lattice inter- 
actions and energetics. Hence, we have also performed 
total energy calculations from density functional the- 
ory (DFT). The calculations were performed within the 
plane-wave implementation iTsj of the local-density ap- 
proximation (LDA) to DFT. We used Vanderbilt ultra- 
soft pseudopotentials|0. We relax only the hydrogen 
molecules inside the primitive cell of the M0F5 struc- 
ture, which contains 106 atoms fl^. A cutoff energy of 
340 eV was found to be enough for total energies to con- 
verge within 0.5 meV/atom. 

The energies of the four adsorption sites are summa- 
rized in Table 1 for two different orientation of the hy- 
drogen molecule. The binding energies are in good agree- 
ment with the experimental finding the " Cup-site" is the 
most energetically stable, followed by the "ZnOa-site". 
The calculated binding energies for the "Hex-site" and 
"Zn02-site" are quite close each other, in agreement with 



TABLE I: The calculated binding energies for four-adsorption 
sites when H2 molecule is parallel and perpendicular to the 
3-foId axis near the adsorption sites. 



the equal population of these sites observed experimen- 
tally. Finally, we also point out that the "Hex-site" and 
"Zn02-site" are further stabilized by the intermolecular 
interactions amongst the adsorbed deuterium molecules. 
For example, each "ZnOa-site" is surrounded by three 
"Zn02-site" due to local 3-fold symmetry. Upon initial 
hydrogen loading, the molecules first go to the biggest 
cavities, namely the "cup-site" and the corner sites (i.e. 
"ZnOs-site"). After that, each hydrogen molecules is 
surrounded by three hydrogen molecules (of the "Zn02- 
site" and "Hex-site"). In this way, the packing of the 
hydrogen molecules are optimized due to both H2-ZnO- 
cluster interaction and the H2-H2 interactions. Finally, 
we note that the adsorption energy of the first three sites 
shows significant anisotropy (about 30 meV) with re- 
spect to the orientation of the H2 molecule. Therefore 
we expect significant splitting of the ortho-para transi- 
tions for H2 molecules adsorbed at these three sites with 
a more isotropic transition for the hydrogen molecule 
at the "Hex-site". These results seem to be consis- 
tent with the available inelastic neutron data on H2 in 
M0F5|^|17,]. We also observe that the exact location of 
the H2 molecule at the adsorption sites depends on the H2 
orientation within 0.2 A , suggesting that a proper treat- 
ment of translation-rotation coupling of the H2 quan- 
tum dynamics is required to explain the inelastic neutron 
datatla- 

So far we have discussed hydrogen loading up to 26 
D2/4Zn, at which point the adsorption sites discussed 
above are almost fully occupied. Figure 5 shows the 
neutron powder diffraction patterns and Rietveld struc- 
tural refinements at two more deuterium loadings of 34 
and 46 D2/4Zn, indicating that M0F5 structure is capa- 
ble of adsorbing more hydrogens. The difference Fourier 
analysis indicated that at these high-coverages, hydro- 
gen molecules form quite interesting nano-cages in the 
cubic-cavities of the M0F5 structure as shown in Fig. 6. 
The first two hydrogen positions listed in Fig. 5 gener- 
ate the cubic and almost spherical D2 nano-cages shown 
in Fig. 6(a) and (b), respectively. Due to the hydrogen 
atoms on the organic linker, the cubic nanocage shown in 
Fig. 6(a) is slightly bent along the edges. For 34 D2/4Zn 
loading, these two cage structures have about equal pop- 
ulation; indicating some disorder. However, with increas- 
ing the hydrogen loading to 46 D2 /4Zn, we observed that 
the cubic cage is destabilized with respect to the more 
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FIG. 5: The neutron powder diffraction patterns (dots), 
Rietveld refinements (solid) and the difference plots (noisy 
line) for high-concentration hydrogen loading, n=34 and n=46 
D2 /4Zn, respectively. In addition to the four adsorption sites 
shown in Fig. 4, additional hydrogen sites, occupancies and 
thermal factors are also given. 






symmetric and exotic looking cage shown in Fig. 6(b). 
The intermolecular distances in these nano-cages are on 
the order of 3.0 A, much shorter than those found in solid 
H2^3;EI- At the maximum coverage of 46D2/4Zn, we 
also determined three additional hydrogen sites which 
are listed in Fig. 5. These hydrogens basically sit on 
the top of the square faces of the nano-cage shown in 
Fig. 6(b), creating quite remarkable 3D interlinked nano- 
cage structures. These results suggest that the MOF host 
lattices may be used as a template to build new artifi- 
cial hydrogen nano-structures, which could have quite 
interesting properties due to the quantum nature of the 
molecules, confinement effects, and small intermolecular 
distances. The structure of solid hydrogen under very 
high-pressures has been a focus of intense research for 
a long time due to theoretical predictions for metallic 
behaviorT^. Hence, we hope that our initial results for 
the high-concentration hydrogen loading will give a dif- 
ferent perspective and direction to this important field of 
research [13 • We note that due to the construction of 
the vanadium sample holder used in this study, we were 
not able to apply high pressures to load the M0F5 sam- 
ple with even more hydrogen molecules. The fractional 
occupancies listed in Fig. 5 indicates that it should be 
possible to insert more hydrogen molecules into the struc- 
ture. We hope to perform such high-pressure studies in 
the near future. 

In conclusion, using Rietveld structural refinement 
of neutron powder diffraction data in conjunction with 
difference-Fourier analysis and first-principles calcula- 
tions, we have determined the hydrogen adsorption sites 
and binding energies in M0F5. We have also discov- 
ered that hydrogen molecules form unique 3D interlinked 



FIG. 6: Two different hydrogen nano-cages (blue and green 
spheres) obtained from Rietveld refinements at 34 and 46 
D2/4Zn loading, respectively. For clarity, only some of the 
ZnOa (red-pink) and CH bonds (grey-black) are also shown. 



nano-cages at high-concentrations of hydrogen loading. 
Surprisingly, we find that MOF host lattice has enough 
space to hold hydrogen molecules up to 10-wt% at low 
temperatures. This implies that by using different or- 
ganic linkers, which make hydrogen desorption difficult 
by narrowing the channels connecting the network of 
nano-pores in MOF, one may be able to engineer these 
materials for practical hydrogen storage at ambient con- 
ditions. These results not only hold the key to optimiz- 
ing MOF materials for hydrogen storage applications but 
also suggest that MOFs can be used as templates to cre- 
ate artificial interlinked hydrogen nano-cages with novel 
properties. 

Acknowledgments: We thank B. H. Toby and Q. 
Huang for assistance with BT-1 and J. Eckert, J. Rowsell, 
and D. A. Neumann for useful discussions. This work was 
partially supported by DOE within the Center of Excel- 
lence on Carbon-based Hydrogen Storage Materials. 



[1] 

[2] 
[3] 
[4] 
[5] 



See the special issue Towards a Hydrogen Economy, by 

R. Coontz and B. Hanson, Science 305, 957 (2004). 

A. Zuttel, Materials Today 6, 24 (2003). 

M. Eddaoudi et. al., Science 295, 469 (2002). 

O. M. Yaghi et. ai, Nature 423, 705 (2003). 

N. L. Rosi, et. al, Science 300, 1127 (2003); please note 

that the hydrogen isotherm reported in this paper is not 



5 



correct; see J. L. C. Rowsell, A. R. Millward, K. Sung 
Park, and O. M. Yaghi, J. Am. Chem. Soc. 126, 5666 
(2004). 

[6] X. Zhao et. al, Science 306, 1012 (2004). 

[7] N. W. Ockwing, O. D. Friedrichs, M. O'Keeffe, and O. 

M. Yaghi, Acc. Chem.Res. 38, 176 (2005). 
[8] Q. Yang and C. Zhong, J. Phys. Chem. B Lett. 109, 

11862 (2005). 
[9] R. Kitaura et. al, Science 298, 2358 (2002). 
[10] J. van Kranendonk, Solid Hydrogen (Plenum, New York, 

1983). 

[11] R. J. Hemley and H. K. Mao, Phys. Rev. Lett. 61, 857 
(1988). 



[12] D. E. Ramaker, L. Kumar, and F. E. Harris, Phys. Rev. 
Lett. 34, 812 (1975). 

[13] The details of sample synthesis, the Rietveld refinements, 
powder patterns, etc can be found at 
littp: //www. ncnr.nist.gov/staff/taner/h2] 

[14] The unit formula for M0F5 is Zn40i3-(C8H4)3. The con- 
ventional unit cell (Fm3m) has eight of these formulas. 

[15] M. C. Payne et. al, Rev. Mod. Phys. 64, 1045 (1992). 

[16] D. Vanderbilt, Phys. Rev. B 41, 7892 (1990). 

[17] M. R. Hartman, T. Yildirim, T. J. Udovic, C. M. Brown, 
to be published (2005). 



